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A NEW METHOD FOR OBTAINING APPROXIMATE SOLUTIONS OF 
THE HYPERBOLIC KEPLER’S EQUATION 

MARTIN AVENDANO, VERONICA MARTIN-MOLINA, AND JORGE ORTIGAS-GALINDO 


Abstract. We provide an approximate zero S(g,L) for the hyperbolic Kepler’s equation 
S — g arcsinh(S') — L = 0 for g G (0,1) and L G [0, oo). We prove, by using Smale’s a-theory, 
that Newton’s method starting at our approximate zero produces a sequence that converges 
to the actual solution S(g,L) at quadratic speed, i.e. if S n is the value obtained after n 
iterations, then \S„ — S| < 0.5 2 _1 |S — S |. The approximate zero S(g,L) is a piecewise- 
defined function involving several linear expressions and one with cubic and square roots. In 
bounded regions of (0,1) x [0, oo) that exclude a small neighborhood of g = 1, L = 0, we also 
provide a method to construct simpler starters involving only constants. 


1. Introduction 

Keplerian orbits can be described by two parameters, the semilatus rectum p > 0 and the 
eccentricity e > 0, from which all other orbital constants can be derived. According to 
Kepler’s laws, the coordinates (x, y ) at time t of a satellite in the standard reference frame 
(central body at the origin and x-axis pointing to the periapsis) are given by 


Elliptic (e < 1) 

Parabolic (e = 1) 

Hyperbolic (e > 1) 

x - j2^(cos(E) e) 

* = 5(i -d 2 ) 

x - I ^(cosh(R) e) 

y vh?^ E ) 

y = pD 

y ~ yjL j smh(H) 

E — esin(E) = M 

D + SSL — M 

esinh(R) -H = M 

to) 

O 

■+o 

1 

■+o 

iJfK 

II 

M - to) 


where y is the gravitational parameter of the central body and to is the time of passage through 
periapsis [8]. The value M, called mean anomaly, is a linear function of the time. E, D and H 
are the elliptic, parabolic and hyperbolic eccentric anomalies respectively, which are functions 
of the time which are related to M through Kepler’s equation (third line of formulas in the 
table above). The eccentric anomalies are needed to compute the coordinates x and y , hence 
the need for a method to solve Kepler’s equation. 

In the case of circular orbits (e = 0), the equation for the eccentric anomaly reduces to E = M, 
and in the case of parabolic trajectories (e = 1), the cubic equation for D can be solved exactly: 

D _ 3 1 3 M + V9 M 2 + 4 + aj 3 M - VdAf- 2 + 4 

Therefore, the real problem resides on the cases 0 < e < 1 and e > 1, where the equation 
involves a combination of algebraic and trascendental functions. 

For simplicity, the hyperbolic Kepler’s equation is usually written as sinh (H) — gH = L, where 
g = - e G (0,1) and L = ^ G (—oo, oo). Since the left side of the equation is an odd function of 
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H , it is enough to consider L E [0, oo). Moreover, assuming that H can be found, the formula 
for the y coordinate requires sinh (H), so it makes more sense to directly find S = sinh (H) 
solving 


( 1 . 1 ) 


fg,L(S) = S — g arcsinh(5) — L = 0. 


Once that is done, we can compute y as y = —-S and, since cosh (H) = \J 1 + sinh 2 (77), x 

as x = 1 + S 2 - e). 

While there are many articles discussing solutions for the elliptic case (see m, m, 0 , 0, m , 
m, uni, among others), the hyperbolic case has received less attention. Prussing [11] and 
Serafin [12] gave upper bounds for the actual solution of the hyperbolic Kepler’s equation, 
which can be used as starters for Newton’s method since f" L > 0. Gooding and Odell [6| 
solved the hyperbolic Kepler’s equation by using Newton’s method starting from a well-tuned 
formula depending on the parameters g and L. Their approach gives a relative accuracy of 
10~ 20 with only two iterations. Although their starter is a single formula that works on the 
entire region (0,1) x [0, oo), it is too complicated to provide an efficient implementation. 

In contrast, Fukushima [5] focused on the efficiency and simplicity of the starter rather than 
trying to find a universal formula, and produced a starter that is defined by different formulas, 
each valid in a stripe-like region. He showed that his starter converges under Newton’s method, 
but not necessarily at quadratic speed. 


Our approach to solve Eq. (1.1) is to use Newton’s method starting from a value S(g,L) that 
is close enough to the actual solution S(g, L ) to guarantee quadratic convergence speed, i.e. if 
S n denotes the value obtained after n iterations, then | S n — S\ < 0.5 2 ’ 1 ~ 1 |S' — S\. We use a 
simple criterion, Smale’s a-test [2] with the constant ao improved by Wang and Han |lj], to 


decide whether a starter gives the claimed convergence rate. Values that satisfy the test are 
called approximate zeros. 

Definition 1.1 (Smale’s a-test). Let / : (a, b) C M a R be an infinitely differentiable 
function. The value z € (a, b ) is an approximate zero of / if it satisfies the following condition 


«(/, z) = P(f, z) ■ 7 (/, z) < a 0 , 


where 

P(f,z) = 

and a 0 = 3 - 2^2 « 0.1715728. 


f(z) 


f'(z) 


7(/, z) = sup 

k> 2 


/ (fc) (; 


k\f'{z) 


1 

k -1 


In an earlier paper [T] , we obtained the following simple approximate zero E(e,M ) for the 
elliptic case. 


Theorem 1.2 ([1]). The starter 


E(e,M) = < 


M 

27r /3 

V 2 

M 

1—e 


if e < 1 /2 or M > 2?r /3 

if e > V 2 and < M < 2ir /3 

if e > l /2 and 7/7 < M < */a 

if e> 1 / 2 , M < and M < 


sj 6 Me 2 


2(1—e) 

vTmP 


otherwise 


is an approximate zero of E — e sin(K) — M for all e G [0,1) and M G [0, 7 r]. 
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Our main result in this paper, proven in Section [2j is a starter for the hyperbolic Kepler’s 
equation that is a piecewise-defined function in eight stripe-like regions. In seven of them, we 
use only linear expressions in g and L, while in the last one we need a cubic and a square root. 


Theorem 1.3. The starter 


S(g,L) = { 


L + 2.30 g 

if 4 — 1.9 g < L 

L + 1.90 g 

if 2.74 - 1.56 g < L < 4 — 1.9 g 

L + 1.56 <7 

if 2.01 - 1.33 g <L< 2.74 - 1.56 g 

L + 1.33 g 

if 1.60 - 1.16 g <L< 2.01 - 1.33 g 

L + 1.16 g 

if 1.32 - 1.02 g <L< 1.60 - 1.16 g 

L + 1.02 g 

if 1.12 - 0.91 g <L< 1.32 - 1.02 g 

L + 0.91 g 

if 1 - \g < L < 1.12 — 0.91 g 


the exact solution of (1 — g)S + g ^ = L if 0 < L < 1 — | <7 
is an approximate zero of f gt L for all g 6 (0,1) and L 6 [0, oo). 



Figure 1. The regions of Theorem 1.3 are represented in different colors. The 
uppermost region (in dark red) is actually unbounded 


Remark 1.4. The exact solution S of the cubic equation (1 — g)S + g 


sf 

6 


L is given by 


S(g,L) 




8(1 -g) 3 

n 3 


+ 



which can also be written as 


S(g, L) = A 


2(1 -9) 

gA 


where A 




8(1 -g) 3 

n 3 



8(1 -g) 3 

n 3 


For any bounded region of (0,1) x [0,oo) that excludes a neighborhood of g = 1,L = 0, 
we give in Section [3] an explicit construction of another piecewise-defined starter using only 
constants. This is important for building look-up table solutions of the hyperbolic Kepler’s 
equation. Finally, we show that such a construction is impossible if a neighborhood of the 
corner is included. 
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Theorem 1.5. For any 0 < e < \ and L max > e' = v ^ a ° £ 1 2 . there is a piecewise constant 

(l-ep 

function S defined in R = (0,1) x [0, L max ] \ (1 — e, 1 ) x [0, e'] that is an approximate zero of 

fg,L- 

Theorem 1.6. Let U C (0,1) x [0, oo) be an open set such that U D {1} x [0,e] for some 
e > 0. For any constant starter So > 0, there exists a point (go,Lo) £ U such that So is not 
an approximate zero of f go j jQ . 


2 . Starters for the hyperbolic Kepler’s equation 


The aim of this section is to prove Theorem 1.3 To do that, we first prove some necessary 
technical results. Then we show in Theorem 2.4 that the solution of the cubic equation 
(1 — g)S + 9_j£_ = L is an approximate zero when 0 < L < 1 — | g. Finally, we study 
in Theorem 


2.5 


the family of starters S = L + ag , and in Corollary 2.6 those with a E 


{0.91,1.02,1.16,1.33,1.56,1.90, 2.30}. 


Lemma 2.1. For all k > 1 , we have 

dfi , 

—r(arcsinh(x)) = (1 + x 2 )^ ~ k Pk(x), 
dx K 

where Pk(x) are polynomials of degree k — 1. The monomials of Pk(x) have all even or odd 
degree: 

Pk(x) = aj^iz*" 1 + a { k %x k ~ 3 + a ( k l 5 x k - 5 + ■■■ 

The leading coefficient is (—1 ) fc-1 (A; — 1)!, the coefficients alternate signs, and the sum of 
the absolute value of all the coefficients is ||Pfc(x)||i = 1 • 3 • 5 • • • (2k — 3) = (2k — 3)!! The 
independent term satisfies 


^( 0 )| 


((k — 2)!!) 2 if k is odd, 

0 if k is even. 


Proof. We proceed by induction. Note that arcsinh^x) = (1 + x 2 )~z, so the Lemma holds 
for k = 1 by setting P\(x) = 1. Assume now that the Lemma is true for some k > 1. 
Differentiating the fc-th derivative, we get 

d tk+1 arcsinh(x) = (1 + x 2 )2~ ( ' k+1 ' > [(1 — 2 k)xPk(x) + (1 + x 2 )P(.(x)\ , 

so we define Pk+i(x) = (l — 2k)xPk(x) + (l+x 2 )P k (x). This shows that Pk+i(x) is a polynomial 
of degree at most k. Moreover, if Pk(x) has only terms of odd (or even) degree, then Pk+\(x) 
has only terms of even (or odd) degree, respectively. A more careful study of the leading 
coefficient of Pk+\ (x) from the recurrence gives 

af +1 ) = (1 - 2k) a™, + (k - l )^ k \ = (-^afl, = (-k)(-l) k -\k - 1)! = (-1 ) k k\, 


showing that the degree of Pk+i(x) is k and that the leading coefficient of Pk+i(x) is (— l) k k\, 
as we needed for the inductive step. The other coefficients of Pk+ i(x) can be also found by 
using the recurrence: 


a 


(fc+i) 

k—2r 


(~ k ~ 2r)4 fc 2 2r _ 1 + (k -2 r + 1 )a^ 2r+1 


rkl (k) (k) 

for r = 1,..., . Since —k — 2r < 0, k — 2r + 1 > 0, and the signs of a k _ 2r _ 1 and a k _ 2r+1 are 

different, we have sgn(a^l l 2 ^) = sgn(a^ 2r+1 ), which alternate by the inductive hypothesis. 
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Finally, for a polynomial with coefficients of alternating signs and all even (or odd) exponents, 
we have ||P fc+ i(x)||i = |Pfc+i(i)l = |(1 - 2k)iP k (i)\ = (2k - l)||P fc (P)||i = (2k - 1 )!! 

According to [7], the power series expansion of arcsinh(x) at x = 0 is 

3 inh(x) = V (-!) n ( 2 n ~ 1 )H 2n+1 = ^ k 

v ^ (2n + l)( 2 n)!! k\ 


arcsmi: 


n> o 

so P 2 n (0) = 0, and P 2 n+ 1 (0) = = (~l) n (2n - l)!! 2 . 

Lemma 2.2. Given a real number i/O and an integer n> 2, then 


k> o 


sup 

k>n 


X 

fc-i 

X 


— 

= max < 1 , 

— 


k 

l 

n 

J 


□ 


Proof. Assume first that |x| < n. Then max{l, I-I"- 1 } = 1 . Moreover, 



X 

1 

k — 1 

X 

sup 

k 

< sup 

— 

k>n 

k>n 

n 


i 

fc-i 


< 1 


and 


Therefore, 


sup 

k>n 


X 

k 1 > lim 

X 

k 

k —^00 

k 


sup 

k>n 


l 

fc-i 


TTTTT lim log Pl-log |fc| 

— gfc-Aoo k — 1 — — p 


= 1 = max < 1 , — 
I n 


Now consider the case Ixl > n. Then max < 1 


-it — Moreover, 

71 . 


sup 

k>n 


k 


k — 1 ^ I I n — 1 

~ I n I 


since the supremum is greater than or equal to the first term, and 



X 

1 

k-1 . 

X 

1 

k-1 . 

X 

sup 

k>n 

k 

< sup 

k>n 

n 

< 

n 


i 

71 — 1 


X 

k-1 

X 

71 — 1 


X 

71-1 

k 


n 

= max < 

n 



since the sequence |^| n_1 is decreasing. Therefore, 

i 

x / 

sup 

k>n 

which concludes the proof. 

The previous two lemmas give us the following upper bound for 7 (f g ,L, S). 
Lemma 2.3. For all n > 2, and (g, L ) e (0,1) x [0, oo), 


□ 


7 (fg,L,S)< 


where 


1 


1 + S 2 


max < a 2 ,..., a n _ 1 , 2 iriaxjl, |5|} max < 1 , 


9 


,n(V 1 + S 2 — g) t 


^k — 


( g\Pk(S)\ 
\k\(Vl + S 2 -g), 


1 

k-1 
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Proof. By Lemma 2.1 


7 (fg,L,S) = sup 
k> 2 


= sup 

k> 2 


fl%s) 


k'-f g , L (S) 

gPk&il + 'S 2 ) 1 -* 


1 


= sup 
k> 2 

5 P fe (5)(l + 5 2 )i- fc 

4-^) 


1 

k -1 


fc! (Vl + S 2 -g 


i 

fe-i 


1 


-— sup 

1 + S 2 k> 2 


^(5) 


A:! ( VI+ S 2 ~g 


1 

fc-1 


1 


1 + S' 2 


max < a 2 ,..., a n _i, sup a*, 

fc>n 


Again by Lemma |2.1 


gmaxjl, |S'|} fc 1 • (2k — 3)!! 

k-1 

g ■ (2k — 2)!! 

k\ (Vl + S 2 — g^j 


k\ (Vl + S 2 -g^ 


i 

fc-i 


= max{l, |5|} 


9 -2 


fc-i 


HVl + S 2 -9 


i 

&-i 


= 2max{l, |5|} 


k V1 + S 2 - g 


i 

fc-i 


By Lemma 2.2 


sup a k < 2max{l, |5|}sup 

fc>n k>n 


9 


= 2 max{ 1, | S\ } max < 


k(Vl + S 2 -g 
9 


i 

k -1 


1 

n —1 


1, 


n ( V 1 + S’ 2 — g 


. 


□ 

Theorem 2.4. The exact solution S of the cubic equation (1 -g)S+g^r- = L is an approximate 
zero of fg t L in the region R\, where 


Ri = |0 < g < 1, 0 < L < 1 — -g 

Proof. For any point (g, L) £ R.\ , the solution S(g, L) of the cubic equation satisfies S £ [0,1]. 
Therefore, 

\f g ,L{S)\= S — g arcsinh(S) — L = ((1 - g)S + |s 3 - l)j - g ^arcsinh(S’) - S + ^S 3 ^j 

= g ^arcsinh(S') — S + ^P 3 ^ < arcsinh(S') — S + ^S 3 . 

Using the previous inequality in the definition of /3(/ 9 i l, S) and the fact that g < 1, we obtain 


P(fg,L,S) = 


f(S ) 


ns) 


< 


arcsinh(S) - S + |S 3 _ (arcsinh(S) - 5 + |S 3 ) \/l + S 2 


1 - 


\A + S 2 


Vl + S 2 -1 
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On the other hand, it follows from Lemma 


2.3 


7 (fg,L,S)< 


1 


s 


1 + 5 2 
1 


max 


max < 2 , 


2(\/T-i-S' 2 — 1) 

2 


for n = 3, g < 1 and |5| < 1 that 

. , 1 \ 

, 2 max < 1 , 


,3(^1+ 5 2 - 1), 


1 + S 2 1 ’ V3(Vl+¥-l)^ I ’ 


since 


2(v / l+5 2 -l) 


< 


-1 


Multiplying the inequalities for /3(/ 9 i l, 5) and 7 (/ s ,l, 5) together, we have that the ce-test is 
satisfied if 


arcsinh(S’) — 5 + |5 3 ^ 

Vi + s 2 (Vi + 5 2 - 1 ) 1 ’VsiVi + s 2 - 1)2 


max < 2, 


< a 0 - 


A standard analytic study of the previous one-variable function shows that it is bounded above 
by ao in the interval 5 E [0,1]. □ 

Theorem 2.5. The starter S(g, L) = L + ag is an approximate zero of f g j J in the stripe 

R 2 ^a) — (0 < <7 < 1, 5 m ; n (ci) 7 L T ag 7 5 max (ci)} 

if Smin(o), 5 max (a) chosen such that 

\a — arcsinh(x)| r , /-„ N ^ 

. — . -max < 1 , x, 2 x(v 1 + x 2 — 1 ) > < cto 

t/TT^VTT^ 2 -i ) 2 l > 0 

/or all x E [5 min (a) ) *S , max(n)] • 

Proof. By definition, we have 


P{fg,L,S) = 


fg,L(S) 




g\a — arcsinh(5)| \a — arcsinh(5)| 


By Lemma 2.3 with n = 2, 


7 (fg,L,S)< 


l + S 2 

Therefore, we obtain that 

<x(f g>L ,S) = P(f g ,L,S)-'y(f g , L ,S) 
21 a — arcsinh(5)| 


1 


max{l, 5} max < 1, 


1 - 


7 / 1 +s 2 


2(v/i + 5 2 -i) 


< 


V / l + 5 2 (v / l + 5 2 -l) 

| a — arcsinh(5)| 


maxjl, 5} max < 1, 


2(v/l + 5 2 -l) 


V / l + 5 2 (v / l + 5 2 -l) 2 
which is less than ao in R- 2 (a) when 5 E [5 m i n (a), 5 max (a)]. 


max < 1, 5, 25( V 1 + 5 2 — 1) > , 


□ 
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Corollary 2.6. The starter S{g, L) = L + ag is an approximate zero of f g x in the stripe 
R“2 (®) — {0 ^ 9 ^ 1) S m i n (u) 7 L -l~ 0,9 ^ S max (o)} 


for the values of a, S m i n (a), S max (a) given in the following table: 


a 

Smin(®) 

Smax(®) 

0.91 

0.99 

1.12 

1.02 

1.12 

1.32 

1.16 

1.32 

1.60 

1.33 

1.59 

2.01 

1.56 

2.00 

2.74 

1.90 

2.73 

4.00 

2.30 

4.00 

oo 


2.4 


S(g, L ) is an app roximate zero of f g .r when 0 < L < 
It is easy to verify that for each 


2.6 


Proof of Theorem m By Theorem 
1 — | g. For the remaining cases, we use Corollary 
a € {0.91,1.02,1.16,1.33,1.56,1.90, 2.30}, the region where S is defined as L + ag is contained 
in i? 2 (a). □ 


3. An analytical study of piecewise constant starters 


Theorem 3.1. The constant starter S(g,L ) = So > 0 is an approximate zero of f g _i in the 
stripe 

S 3 (S 0 ) = {o < g < 1, So - g arcsinh(S 0 ) - A(S 0 ) < L < S 0 - g arcsinh(S 0 ) + A(S 0 )} , 
where 


A(S 0 ) = 


ifO<S 0 <^, 


V3aoS$y/l+Sg 

2(V 1 + S 0+ 1 ) 3/ ""’ 

aggp minjl.Soj-y/l+gg X C V7 

2(Vl+ S 0+ 1 ) ’ 1 ° _ 3 ' 


The constant starter S(g,L) = 0 is an approximate zero of f g L in the region 

\/3a 0 (l -j)H I 


Rs(0) = <0 < g <1,0 < L < min { a 0 (l - g ), 


1 

g 2 


Proof. Consider first the case So > 0. By Definition 0 and the fact that g < 1, 

P{fg,L,S) = 


f 9 ,L(S) 


So — garcsinh(So) — L 

f’ 9 ,L(S) 


1 9 

yj l + s a 


< 


| So — 9 arcsinh(So) — L\ 1 


1 - 


V /l+s o 


= ^2 I'S'o - S'arcsinh(So) - L\ \Jl + S'l ( y^l + S^ + 1 j . 


On the other hand, it follows from Lemma |2.3| with n = 3 that 

So 


7 (fg,L,S)< 


1 

TTs^ 

2 

l + 5 n 2 


max 


2(v / TTS 2 - 1) 
max{l, Sq} max < 1 


, 2 max{l, Sq} max < 1, 


3(^1 + Sq - 1) 


3{\/l + Sq — 1) 
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We will now distinguish between So > 1 and 0 < So < 1. In the first case, 7 {f g ,Li S) < 37^1 
so the a-test follows from 

2 ^ sj\ + Sq + l) |S 0 - 5 arcsinh(S 0 ) - L\ 


Sov/l + Sg 


< ot 0 , 


which is equivalent to 


IS 0 - 5 arcsinh(So) - L\ < ^ 1 + ^ = A(S 0 ). 

2 (V 1 + S 0 + *) 

When 0 < Sq < 1, 7 (f g ,L, S) < y+S 2 max | 1) f so a ~test is satisfied if 


| So — ^ arcsinh(So) — L 
which is equivalent to 


2 (v 1 + 7 +1) 


sjv'iTs? 


max < 1 , 


3(^/1 + Sq - 1) 


< CKQ, 


|Sq — g arcsinh(So) — L\ < min < 


aoS 0 2 v / TTs2 V3aoS 0 3 yrTs2 


2 (v /l + 5 o + 1 ) 2 ( Wl + S2 + 1 


= 


\/3aoSg-y/l+5o ,p „ 

2 (V 1 + S 0+ 1 ) 3/2 


<S 0 < ^ 


aoSi^T+Sj if ^ < S 0 < 1 

2(71^52+1) 3 - 


= A(S 0 ). 


Therefore, the a-test is satisfied in the whole region i? 3 (So). 
When S(< 7 , L) = So = 0, we have 

fg,L(S) L 


P(fg,L,S) = 


fU S ) 


1 ~9 


and by Lemmas 2.1 and |2.2| 


7 ifg,L,S) = sup 
k> 2 


f (k) { 0) 


< sup 
k> 3 


k\(l-g) 
9 


= sup 

k> 2 


gPk{ 0 ) 


k\(l-g) 


= sup 

k> 3 
odd 


9 ((k — 2)!!) 2 


k\{l- g) 


k{ 1 - 5 ) 

Therefore, the a-test is satisfied if 

L 


= max < 1 , 


3(1-5)/' 


1 ~9 

which is equivalent to ( 5 , L) G i? 3 ( 0 ). 


max < 1 , 


3(1-5) 


< ot 0 , 


Remark 3.2. The region i? 3 (0) of Theorem 3.1 can be written as 

i? 3 (0) = <fo < 5 < ^,0 < L < a 0 (l -5) j u \ 7 < 5 < 1,0 < L < 


\/3a 0 (l - 5)2 


□ 
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Figure 2. The region Rs(0) of Theorem 3.1 and Remark 3.2 is represented in red 


Lemma 3.3. The function A (So) of Theorem 3.1 is strictly increasing for So > 0. 
Proof. When 0 < Sq < ^, then 


A(S 0 ) = 


v^QpSoy/1 + Sq _ VSao e f 
2(^1 + ^ o 2 + !) 3/2 ” 


S n 2 1 - 


1 T *Sn 


1 - 


"\/l + Sq + 1 


When So > then 


A / C s a 0 S 0 minjl, S 0 } V 1 + Sq «o c n . 1 \ 

A(So> - KvrrwTi) = mm{1 ' s »> ' TttWTT j ■ 

Since both expressions are the product of constant and strictly increasing functions and A (So) 
is continuous at So = y, then A(Sq) is strictly increasing. □ 


Proof of Theorem 1.5 We decompose the region of the theorem as 


(0,1 - e) x [0, s'] U (0,1) x [s', L max ] 


B 


From Remark 3.2, we have that A C i?3(0). 

Define the sequence So = s' and S* + i = Sj + 2A(Sj), for i > 0. Since A is a positive and 


strictly increasing function by Lemma 3.3 the sequence Sj is strictly increasing and satisfies 
Si > So + 2iA(So) —> oo as * —> oo. 

Since lim (x — arcsinh(x)) = oo, there exists a large enough N for which Sn — arcsinh(Sjv) > 

#—>•00 

N 

L m M . Once we prove that B C C = |^J R^Sf), we will have that A U B is covered by the 

i =o 

regions where the starters 0, So,..., Sn are approximate zeros of / 5 y In particular, this will 
provide the piecewise constant starter S that we looked for. 

It only remains to prove that B C C. Note first that i? 3 (Sj) are stripes whose union is 


C = (0 < g < 1 , Sq — A(Sq) — g arcsinh(So) < L < Sn + A(Sn) ~ 9 arcsinh(STv)} 


since the lower boundary of R-j,{Si + \) is a segment that is always below the upper boundary of 
Rs(Si). Indeed, it is enough to see that the endpoints of the first segment, (0, Sj+i — A(S i+1 )) 
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and ( 1 , Sj+i — A(5j + i) — arcsinh(5i + i)), are below the endpoints of the second segment, ( 0 , Si + 
A (Si)) and (l,5j + A (Sf) — arcsinh(S'j)), respectively: 

S i+ 1 - A(S i+1 ) <Si + A (Si) <*=► A (Si) < A (S i+1 ), 

Si + 1 — A(5j + i) — arcsinh(S l j_ l _i) < Si + A (Sf) — arcsinh(S'j) <^=i> 

A (Si) + arcsinh(5j) < A(5j+i) + arcsinh(Si+i), 

which is true because A and arcsinh are both strictly increasing. 

On the other hand, the vertices of the rectangle B are in C, which means that B C C. 
Indeed, it follows from arcsinh(e / ), arcsinh (Sat), A(e'), A(Sjv) > 0 ; e' < L max by hypothesis 
and Sjy — arcsinh(S'jv) > L max by construction, that 

s' — A(e / )— arcsinh(e / ) < e' — A(e / ) < e' < L max < 

< S tv — arcsinh(S’jv) < Sjv + A(Sj\r) — arcsinh(S^r) < $N + A(SV)- 

This implies 

s' - A(e') < e' < S N + A (S N ) => (0, e') G C, 

e' — A(s: / ) — arcsinh(e / ) < e' < Sn + A(Sjv) — arcsinh(S'jv) => (l,^) G C , 

J - a( £ ') < •^max < S'at + A(Sjv) =>• (0, ^max) ^ C*? 

s' - A (s') - arcsinh(e') < L max < Sat + A(Sjv) - arcsinh(SW) => ( 1 , L max ) G C, 
concluding the proof. □ 


Proof of Theorem M We proceed by contradiction, i.e. we assume that Sq is an approximate 
zero of fg t L for all (g, L ) G U. By definition, this means that 


(3.1) 


So — g arcsinh(So) — L 


gPk(So) 

x/l + (V 1 + Sq ~ g) 

k> 2 

k\ (a/ 1 + 5*0 - 5 ) 


< O' o 


for all ( g,L ) G U. We will derive a contradiction from (3.1) for all Sq > 0. 
Case Sq = 0. Inequality (3.1) is equivalent to 


L 

- -sup 

1 ~ 9 k >2 


gPk{ o) 


fc!(i - g) 


< ao- 


Using Lemma 2.1, we obtain 


«0 > 


L 

- -sup 

l — o fc>3 

odd 


g{{k~ 2)!!) 2 

fc-i x, 

g 

fc!(l- ff ) 

> sup 

1 — g k >3 

k(k-l)(l-g) 


1 

fc-1 


> 




V6(l -5)5 

This implies that L < in U. Therefore L < 0 in U D {g = 1} = {1} x [0, el, which 

is impossible. 

For the rest of the cases, we take limit as g —> 1~ and then limit as L —> 0 + , obtaining 


(3.2) 


Sq — arcsinh (Sq) 


V 1 + 'S'o (V 1 + 'S'o ~ *) 


sup 

k> 2 


Pk(S 0 ) 


k\ 


1 + S 2 - 1 


1 

fc -1 


< d 0 . 
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Case 0 < So < 1.598. Inequality (3.2) implies that 

(So - arcsinh(5o))(\/l + + l) 2 ^ 2 (So - arcsinh(So)) 

a o > --— - - > 


25 0 3 yrT5j 


Sq 3 


because (t + l) 2 /t > 4 for all t > 0 . Therefore, h(So) = arcsinh(Sb) — So + -yr 1 > 0 , which is 
false in the interval (0,1.598) because h has only one critical point there, which is a minimum, 
and h(0), h(1.598) < 0. 

Case 1.598 < Sq < 3.1. Inequality (3.2) is equivalent to 

(Sp - arcsinh( 6 , 0 ))(y / l + Sg + 1 )§ ^2S$ - 1 

V6S’ 0 3 v / TTsf 

Since l'V /l+l ^o increasing in this interval, then 


< a 0 . 


ao > 


So - arcsinh(iS’o) (Vl + 1.598 2 + l)i V2 • 1.598 2 - 1 




Therefore, 


So — arcsinh(5’o) 


So 2 


< 


1.598V1 + 1.598 2 
V6a 0 1.598Vl + 1.598 2 


< 0.13, 


(Vl + 1.598 2 + 1)5\/2^T5982TT[ 
and thus So — arcsinh(So) — 0.13Sq < 0 . This is false since the function h(So) = So — 
arcsinh(So) — 0.13Sq has an absolute minimum at So = 3.1 and h( 3.1) > 0 . 

Case 3.1 < So < 9.62. Inequality (3.2) is equivalent to 

4 


(3.3) 


(So - arcsinh(So))(\/l + Sq + l) 3 (6Sq - 9 )3 


< ao- 


V24S*y/T+s$ 

On the one hand, 

( v / ITS i 2 + 1)1 (6Sg-9 )3 > (Vl + 9.62 2 + 1)1(6 • 9.62 2 - 9)1 


> 2.064, 


^oVi + ^o 2 9.62V1 + 9.62 2 

since the expression on the left is decreasing. On the other hand, S ° aiCh ( lll 4- 5 ’o) > 0.24 because 


^0 


h(S 0 ) = So — arcsinh(S'o) — 0.24S 0 3 is increasing in the interval, so h(So) > h( 3.1) > 0. 
Substituting these inequalities in Eq. (3.3), we obtain 

(S 0 - arcsinh(S'o))( v / TT5f + l)t(6Sg - 9)5 2.064-0.24 nif71(7 ^ 

ao > -V- 2 --- > - 57 ^- > 0.1717 > ao, 

V24S 3 v / TTSj V 24 

which is a contradiction. 

Case So > 9.62. Using that > 1 anc [ that Vl + Sq ~ 1 < So i n the interval, we 

_ Vi+'-’o 

obtain from Eq. (3.2) that 


ao > 


1 

Pk(So) 

k-1 

1 

Pk(So) 

/ / - \ 0 U P 

2 (Vl + " l) k ^ 2 

*! (v/i + sg - 1) 

— OO & 

^0 k >2 

k\S 0 


1 

k-1 


This implies that, for all k > 2, 


\Pk(S 0 )\ < (2a 0 ) k ~ 1 k\S k . 
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By Lemma 


2.1 


the leading term of Pk(So) is ±{k — 1)!5 ( 


k—1 


maximum of (2k — 3)!! and it has no term of degree k — 2, so 


its coefficients add up to a 


( 2 a 0 ) k - l k\S% > |Pfc(5 0 )| >{k- 1 )!S * _1 - ( 2 k - 3 )llS k ~ 3 , for all k > 2 , 


which is equivalent to 

(2a 0 ) fc_ 1 A;5'o - 5 0 2 + > 0 , for all k > 2 . 

This implies that 

h k (S 0 ) = {2a 0 ) k ~ 1 kSl - S$ + 2 fc “ 1 > 0, for all k > 2 . 


The degree 3 polynomial h k has two local extrema, a local maximum at So = 0 where h k ( 0) = 
2 fc_1 > 0 and a local minimum at So = r m \ n (k) = ]k _ 1 . Since hk(r m \ n (k)) < 0 if and 

only if k > 5, h k has two positive roots when k > 5, namely n e ft(^) and ?n ght(&)- Therefore, 
hk{So) > 0 is equivalent to So G [0, n e ft(&)] U [bright {k) , °o) when k > 5. 

We will obtain a contradiction by showing that So G (neft(^’)i r right(^)) for some k > 5. This is 
necessarily true since U k > 5 (r Ieft (k), r right (k)) = (r left (5), oo) D (r min (5),oo) D [9.62, oo) 9 S 0 . 
It only remains to be proven that Ufc> 5 (n e ft(P),'/"right(^)) = (Heft (5), oo). This follows from 
Tright(fc) > r min (k) = 3k ( 2 a^h-i —t oo as k -> oo, and the fact that n e{t (k + 1 ) < r rig ht(A:) 
for all k > 5. Indeed, since hy. + \(r m [ n (k)) < 0 and r m i n (fc) < r r \ph t(k), then r’i e ft(fc + 1) < 
r m in(k) < 7-right(&)■ □ 


4. Conclusions 


We provided in Theorem 1.3 a very simple starter S(g,L ) for the hyperbolic Kepler’s equa¬ 


tion which uses a single addition and multiplication when L > 1 — and a few arithmetic 
operations, a cubic and square root otherwise. Our starter can be implemented efficiently on 
modern computers and has guaranteed quadratic convergence speed from the first iteration, 
so the relative error becomes negligible after just a few iterations for any value of g G ( 0 , 1 ) 
and L G [0, oo). 

If an even simpler starter is necessary, we showed in Theorem |1.5| that it is possible to produce 
a piecewise constant starter in any bounded region of ( 0 , 1 ) x [ 0 , oo), excluding a neighborhood 
of the corner g = 1, L = 0. This starter also has quadratic convergence rate and its evaluation 
requires no operations, so it can be efficiently implemented as a look-up table. 
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